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ABSTRACT 

We have calculated the mass accumulation efficiency during helium shell flashes on white 
dwarfs (WDs) of mass 0.7, 0.8, 0.9, 1.0, 1.1, 1.2, 1.3, and 1.35 Mq for the helium accretion rates 
logMne {Mq yr~^) — —7.4 to —5.8. This efficiency is a crucial factor for binary evolutions of 
Type la supernovae. For less massive WDs (< 0.8 Mq) no wind mass loss occurs and all the 
accreted mass accumulates on the WD if the Roche lobe size is large enough. The efficiency 
takes the minimum values in between 1.1 and 1.2 Mq WD for a given mass accretion rate and 
increases in both less and more massive WDs. The mass accumulation efficiency is larger than 
0.5 for logMno > -6.72 in all the WD masses. 

Subject headings: binaries: close — novae — stars: mass loss — supernovae: general — white dwarfs 



1. Introduction 

Recent findings on the acceleration of the uni- 
verse are based on the brightness of Type la su- 
pernovae (e.g., Perlmutter et al. 1999; Riess et al. 
1998). However, the diversity of the brightness has 
not been fully understood yet, which may be re- 
lated to their binary evolution paths (e.g., Umeda 
et al. 1999). In order to elucidate the physics 
of Type la supernovae, every evolutionary path 
to Type la supernovae should be examined, even 
though its production rate is very rare. In these 
paths, the response of white dwarfs (WDs) to hy- 
drogen/helium matter accretion is a key process 
and, especially, the efficiency of mass accumula- 
tion is crucial. 

Recent binary evolution scenarios of Type la 
supernovae include a WD accreting matter lost 
by a companion (Hachisu, Kato, & Nomoto 1996, 



1999a; Hachisu et al. 1999b; Han & Podsiadlowski 
2004; Langer et al. 2000; Li & van den Heuvel 
1977), although the merger model of double degen- 
erates has not been fully rejected yet (e.g., Livio & 
Reiss 2003). In these scenarios, the WD grows in 
mass and explodes as a Type la supernova when it 
reaches the Chandrasekhar mass limit. When the 
accretion rate of hydrogen is much smaller than 
~ 1 X IO^^Mq yr^^, hydrogen shell flash triggers 
a nova outburst. Almost all of the accreted matter 
is blown off or even the WD is eroded by convec- 
tive dredge-up (e.g., Prialnik 1986). When the 
accretion rate is as large as ~ 1 x 10"'' Mq yr~^ 
or more, the shell flash is relatively weak or hy- 
drogen burning is stable so that most part of the 
accreted matter is processed to helium and then 
accumulates on the WD. 

Another possible way to Type la supernovae is 
an helium accretor; the WD accretes helium from 
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companion helium star and grows in mass to the 
Chandrasekhar mass limit (e.g., Yoon & Langer 
2003). In such helium shell flash occurs to 

develop a nova outburst when the mass accretion 
rate is smaller than a critical value, which varies 
from 1 X 10~^ to 1 xlO~^M0 yr~^, depending on 
the WD mass. If the helium shell flashes occur, 
the wind mass loss carries away a part of the en- 
velope matter, which reduces the growth rate of 
the WD. Yoon & Langer (2003) have calculated 
binary evolution that contains a helium accret- 
ing WD. In their calculation, the effect of wind 
mass loss is partly taken into account based on an 
empirical formula of Wolf-Rayet star mass loss, 
whereas Han & Podsiadlowski (2004) have used 
the mass accumulation efficiency of helium accre- 
tion calculated by Kato & Hachisu (1999) for a 
1.3 Mq WD. In this way, the mass accumulation 
efficiency is a crucial factor for evolution of the 
binary. 

Ivanova & Taam (2004) have calculated evo- 
lution of WD binaries in a phase of thermal 
timescale mass-transfer. The crucial factor of the 
binary evolutions is the reaction of WDs against 
high mass transfer rate, that is, the mass ac;cre- 
tion efficiency for hydrogen/helium shell flashes. 
Taam (2004) addressed the mass accumulation 
efficiency of helium novae is one of the three fun- 
damental factors that should be clarified in order 
to advance the evolutionary scenarios to the next 
step. The mass accumulation efficiency for helium 
shell flashes has been calculated but only for a 
1.3 Mq WD (Kato & Hachisu 1999). 

From the observational point of view, such a 
helium shell flash should be observed as a helium 
nova outburst. Recently, in late 2000, V445 Pup- 
pis was discovered at its outburst stage. Ashok 
& Banerjce (2003) suggested, from its infra-red 
spectrum features, that V445 Pup is the first iden- 
tified helium nova. Kato & Hachisu (2003) calcu- 
lated light curves of helium novae for various WD 
masses to fit them with the observational data. 
They concluded that the light curve can be well 
fitted by a helium nova on a WD having its mass 
more massive than 1.33 Mq. This massive WD 
indicates that the WD is now growing in mass 
even for such a violent helium shell flash. Kato & 
Hachisu (2003) suggested that V445 Pup is a pro- 
genitor of Type la supernova. Therefore, we are 
forced to realize the importance of mass loss/mass 



accumulation during helium shell flashes on the 
WD. 

Kato, Saio, & Hachisu (1989) calculated the 
evolution of helium shell flashes on an 1.3 Mq WD. 
Kato & Hachisu (1999) have recalculated the he- 
lium shell flashes with OPAL opacity and obtained 
the mass accumulation efliciency, i.e., the ratio of 
the envelope mass accumulated after one cycle of 
helium shell flash to the ignition mass. These two 
works are done only for a 1.3 Mq WD and no re- 
sults have been presented for other WD masses. 
If the accumulation efficiency is very small for less 
massive WDs, the evolution scenarios should be 
drastically changed. In this Letter, therefore, we 
have calculated mass accumulation efliciency of 
helium shell flashes for various WD masses. 

2. Input Physics 

The decay phase of helium nova outbursts is 
followed based on an optically thick wind theory 
(Kato & Hachisu 1994). The structure of the 
WD envelope is calculated by solving the equa- 
tions of motion, continuity, energy transport, and 
energy conservation. The details of computation 
have been already presented in Kato & Hachisu 
(1999) for helium shell flashes on a 1.3 Mq WD. 

Here, we examine other various WD masses, 
i.e., 0.7, 0.8, 0.9, 1.0, 1.1, 1.2, and 1.35 Mq. The 
radius of the WDs are assumed as given in Ta- 
ble 1 bccaiisc of the following reasons: In time- 
dependent calculation (Kato et al. 1989), the 
radius of the nuclear burning shell moves back 
and forth during one cycle of a helium shell flash. 
In the 1.3 Mq WD case, we assumed the ra- 
dius as a time-averaged mean value of the radius 
(see, e.g., Kato et al. 1989; Kato & Hachisu 
1999). For WD masses other than 1.3 Mq, how- 
ever, no time-dependent calculation of helium shell 
flashes is presented. Therefore, we estimate the ra- 
dius from an geometrical mean of the two radii, 
i.e., the naked C-l-0 core radii of two steadily- 
accreting WD with the mass accretion rate of 
Mue = 1-0 X 10"^ and 1.0 x 10"^ Mq yr-i(private 
communication with H. Saio 2003). 

After the onset of a shell flash, helium is pro- 
cessed to carbon and oxygen. Helium burning pro- 
duces nuclear energy, a large part of which is con- 
sumed for the envelope matter to be push upward 
against the gravity. The nuclear energy release is 
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almost comparable to the gravitational potential, 
and a large part of the envelope is processed in 
very early stage of explosion, i.e., before the ex- 
pansion of the envelope. The processed matter is 
then mixed into the whole envelope by convection. 
Therefore, we assume the composition of the enve- 
lope to be uniform with the values of Y and C + 
as listed in Table 1. Here, we assume X = 0.0, and 
Z = 0.02, where Z includes C and O in the solar 
abundance ratio. The value of C + O = 0.5 for the 
1.3 Mq WD is taken from Kato & Hachisu (1999). 
Changing the ratio of carbon to oxygen with the 
total mass ratio constant (C + O =const.) hardly 
changes the result. OPAL opacity is used. 

3. Mass Accumulation Efficiency 

Figure 1 shows the photospheric temperature 
Tph, the photospheric wind velocity Vjph, the pho- 
tospheric radius i?ph, the wind mass loss rate 
Afwind {dotted) ^ and the total mass decreasing rate 
of the envelope A'/tot = -^'^wind + Afnuc, {solid), here 
the mass decreasing rate owing to nuclear burning 
is Mnuc = Lnuc/{£HcY), whcrc Lnuc is the nuclear 
luminosity. The solution of 1.3 Mq is already pre- 
sented in Kato & Hachisu (1999). At the max- 
imum expansion of the envelope after the onset 
of a shell flash, the star reaches somewhere on the 
curve, depending on the envelope mass at ignition, 
and moves leftward in time. The wind mass loss 
stops at the point marked by the small open cir- 
cles. After that, the star continues to move left- 
ward owing to nuclear burning. Helium nuclear 
burning stops at the end of each curve. In the 
case of 0.8 Mq WD, we have obtained only a very 
short sequence, because of a numerical difficulty 
(see Kato & Hachisu 1994, for more details). No 
winds occur in the 0.7 Mq WD. 

Figure 2 shows the mass accumulation effi- 
ciency r^Ho defined by the ratio of the processed 
matter remaining after one cycle of helium shell 
flash to the ignition mass. The amount of mat- 
ter processed or lost by the wind arc calculated 
from the wind mass loss rate and nuclear burn- 
ing rate in Figure 1 (see Kato et al. 1989, 
for details). Here, we use the relation between 
the ignition mass and the helium accretion rate 
(private communication with H. Saio 2004) ob- 
tained based on a steady state in which non- 
homologous term of gravitational energy genera- 



tion is neglected (Nomoto 1982). For the case of 
0.7 Mq, no wind occurs and all of the envelope 
matter is processed to accumulate on the WD, if 
the binary separation is large enough for the ex- 
panded envelope to reside in the Roche lobe. In 
low accretion rates (logMHe(-^0yr~^) < —7.6) a 
helium detonation may cause a supernova explo- 
sion. 

The value of rjue in Figure 2 is approximated 
by the following algebraic form. 



-0.115(logMHc + 5.7)2 ^ ^ 01, 
?7He=<; (-7.4 < log MHe < -6.05) (1) 

1, (-6.05 < log MHe) 



for 1.35 Mq WD 



-0.175(logMHe + 5.35)2 ^ 03, 
me = S (-7.35 < logMHo < -5.83) (2) 

1, (-5.83 < log MHe) 

for 1.3 Mq WD, 



0.541ogMHe + 4.16 

(-7.06 < logMne < -5.95) 
??He = <! -0.54(logMHc + 5.6)2 _^ ^ 01, (3) 
(-5.95 < logMne < -5.76) 
.1, (-5.76 < logMne) 

for 1.2 Mq WD. This equation is well fitted also 
to 1.1 Mq WD. 



-0.35(logMHc + 5.6)2-^1.01, 
rme=\ (-6.92 < log MHe < -5.93) (4) 
1, (-5.93 < logMne) 

for 1.0 Mq WD, 

( -0.35(logMHc + 5.6)2 _^ i QY, 
me=l (-6.88 < log Mho < -6.05) (5) 
[l, (-6.05 < logMne) 

for 0.9 Mq WD, 



-0.35(logMHc + 6.1)2-^1.02, 
»?He=S (-6.5 < logMne < -6.34) (6) 

1, (-6.34 < logMne) 
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for 0.8 Mq WD. 

When the WD mass is increased, winds begin to 
blow at ~ 0.8 Mq. The winds become stronger as 
the WD mass is further increased. The accumula- 
tion efficiency jjHe then decreases because of much 
stronger wind mass loss. The 1.1 and 1.2 Mq WDs 
take the smallest values of ?7He for all accretion 
rates. For more massive WDs, r/He increases with 
the WD mass because the nuclear burning rates 
are much higher than the wind mass loss rates ow- 
ing to the strong gravity. As a result, most of the 
envelope mass accumulates on these more massive 
WDs. 

4. Discussion and conclusions 

The accumulation efficiency has different values 
if we assume a different chemical composition or a 
WD radius. As the C + O content becomes larger, 
the mass accumulation ratio increases because the 
wind becomes weak but the nuclear burning rate 
increases. For example, if we increase C + O from 
0.5 to 0.6 in the 1.3 Mq case, r/He increases by 
0.08, and from 0.3 to 0.4 in the 1.2 Mq WD, ?7He 
increases by 0.04 — 0.05. If we assume a smaller 
WD radius, the wind becomes stronger, and we 
get a smaller accumulation efficiency TjHe- For ex- 
ample, if we take the Chandrasekhar radius for the 
1.2 Mq WD as the smallest limit, T^He decreases 
by 0.07 at logMne (Mq yr'^) = -6.8 and by 0.12 
at logA^Hc {Mq yr-i) = -6.0. 

For the 0.6 Mq WD, no wind occurs. All the ac- 
creted mass is processed to carbon and oxygen and 
accumulated on the WD. We have not presented 
solutions of the 0.6 Mq WD in Figure 1, because 
the WD radius changes largely depending on the 
mass accretion rate between log A/hc [Mq yr"-'^) = 
—7 and —6. In addition, a density inversion layer 
appears around logT (K) ~ 5.2, which causes a 
numerical difficulty of the computation. 

We have assumed no dredge-up of WD mate- 
rial into the envelope. If WD material under the 
accreted helium layer is mixcid into the envelope, 
the helium content Y is reduced, which increases 
the mass decreasing rate owing to nuclear burn- 
ing becomes larger, because Mnuc fx Y~^. On 
the other hand, the wind mass loss becomes weak 
due to the decrease in opacity, i.e., Mwind becomes 
smaller. Therefore, mixing of WD matter has ef- 
fects to increase the mass accumulation efficiency. 



because 77 ~ Mnuc/ {Mnuc + M^ind) ■ 

For example, in the case of I.3M0, if the same 
amount of WD mass is mixed to the envelope 
mass, i.e., F is ~ 0.5 at ignition. It reduced much 
more in the course of rising phase, because a time 
dependent calculation (Kato, Saio and Hachisu 
1989) shows that Y decreases from 0.98 to 0.5 
in the rising phase. As the gravitational poten- 
tial is comparable to nuclear energy generation 
as shown in table 1, a large part of helium is 
consumed to lift the whole envelope up against 
the gravity at the beginning phase of outburst. 
Therefore, Y is reduced to ~ 0.1 or less. The 
change of Y causes the change in nuclear burn- 
ing rate and the opacity. As a results, the wind 
mass loss does not occur for the envelope less than 
1.5 X 1O~^M0, which correspond to the mass ac- 
cretion rate logM(M0yr~^) = —6.51 i.e., rjcs is 
almost 100 percent at logM > —6.51. Even the 
wind occurs, the mass loss rate is as low as 1/10 of 
nuclear burning rate and almost all of the envelope 
mass will accumulate on the WD. 

It is not known whether WD material is 
dredged up or not, or how much amount is dredged 
up. In the case of hydrogen shell flashes, en- 
richment of WD material is observed in ejecta of 
classical novae, which is an evidence of dredged- 
up of WD material. However, no heavy element 
enrichment is observed in recurrent nova ejecta. 
Therefore, we may conclude that the dredge up 
mechanisms do not effectively work if the recur- 
rence period is short (or the mass accretion rate 
is high). 

In the case of He novae, only one object, V445 
Pup, is observed so far. There arc many bright C 
lines in outburst phase, but we don't know this 
carbon is originated from WD material or pro- 
cessed helium. The WD mass is estimated to be 
> I.33M0, and the recurrence period to be sev- 
eral tens of years (Kato and Hachisu 2003). Such a 
very massive WD itself may suggest that the WD 
is now growing, but not that dredge up mecha- 
nisms are working. 

Cassisi ct al. (1998) have criticized Hachisu et 
al.'s (1996) model in which accreting WDs grow 
in mass and reach the Chandrasekhar mass limit. 
The main point of their criticism is that Hachisu 
et al. (1996) have neglected the interaction be- 
tween the donor and the extended envelope during 
helium shell flashes. Cassisi et al. (1998) specu- 



4 



latcd that most of the envelope mass is lost from 
the binary by the effect of the companion and, as 
a result, the WD cannot grow from ~ IM© to 
the Chandrasckhar mass limit. If so, the Type la 
scenario proposed by Hachisu et al. (1996) is se- 
riously damaged. However, as already discussed 
in Kato & Hachisu (1991a, b, 1994), the viscous 
heating by the companion motion in the envelope 
is not effective when the wind blows. It is be- 
cause the wind velocity is much faster than the 
companion motion and the wind quickly escapes 
from the binary system with almost no interac- 
tion between them. For WDs more massive than 
0.8 Mq, the wind is fast enough and therefore we 
can neglect effects of interaction between the wind 
and the companion's orbital motion. Once the op- 
tically thick wind blows and its velocity is faster 
than the orbital motion, we are able to estimate 
the accumulation efficiency only from the wind so- 
lutions. 

Cassisi et al. (1998) adopted Los Alamos opac- 
ity that shows much smaller values (about one- 
third) than those of OPAL opacity around the 
temperature of logT ~ 5.2. The envelope solu- 
tions for Los Alamos opacity arc largely differ- 
ent from for OPAL opacity when the photosphcric 
temperature is lower than this value (i.e., logTph 
(K)< 5.2). When OPAL opacity is adopted, we 
obtain a more massive envelope mass and a higher 
temperature at the hydrogen burning region for 
the same photosphcric radius. As a result, we 
have a higher nuclear burning rate for a higher 
temperature. A larger part of the envelope mass 
is quickly processed into C + O and accumulates 
on the WD. Therefore, we conclude that a higher 
mass accumulation efficiency is resulted even when 
the strong wind blows. 

Our main results are summarized as follows: 
The mass accumulation efficiency of helium 
shell flashes, TjHe, is a crucial factor for binary 
evolutions of Type la supernovac. For WDs 
more massive than 0.8 Mq, helium shell flashes 
trigger an optically thick wind mass loss. We 
have calculated r/Ho for 0.8 — 1.35 M© WDs; 
more than 50% {rjue > 0.5) is accumulated on 
the WD after one cycle of helium shell flash for 
logMne {Mq yr"i) > -6.72 on the WD masses of 
0.8-1.35 Mq. For WDs less massive than 0.7 Mq, 
no wind mass loss occurs and all the accreted mat- 
ter remains on the WD if we can neglect the effect 



of the Roche lobe overflow. 
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Fig. 1. — Top: Total envelope mass decreasing 
rate, i.e., nuclear burning + wind mass loss {solid 
curve), and wind mass loss rate (dashed curve) in 
units of Mq yr~^, Second: Photospheric radius 
i?ph (cm), Third: wind velocity Vph (cm s~^). Bot- 
tom: temperature Tph (K), against the helium en- 
velope mass, AMue, in units of Mq. The wind 
mass loss does not occur in left side of the point 
marked by a small open circle (an enlarged circle 
for 1.0 Mq) The WD mass is 1.35, 1.3, 1.2, 1.1, 
1.0, 0.9, 0.8, and 0.7 Mq, from left to right. 



This 2-column preprint was prepared with the AAS lAlJ^X 
macros v5.2. 
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Table 1 
Model parameters. 



MwD 

(Mo) 


log RwD 
(«©) 


GMwD/iiwD? 


C + O 


Y 


0.7 


-1.74 


0.12 


0.06 


0.92 


0.8 


-1.84 


0.17 


0.1 


0.88 


0.9 


-1.93 


0.24 


0.13 


0.85 


1.0 


-2.01 


0.32 


0.2 


0.78 


1.1 


-2.09 


0.43 


0.23 


0.75 


1.2 


-2.19 


0.60 


0.3 


0.68 


1.3 


-2.33 


0.88 


0.5 


0.48 


1.35 


-2.44 


1.18 


0.7 


0.28 
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accretion rate dM/dt^^ (Mg/yr) 

Fig. 2. — Mass accumulation efficiency, rynei is 
plotted against helium mass accretion rate for a 
0.7 Mq WD (horizontal solid line), 0.8 and 0.9 Mq 
WDs {thin solid line), 1.0 Mq WD {dash-three- 
dotted line), 1.1 Mq WD {dashed line), 1.2 and 
1.3 Mq WDs {thick solid line), and 1.35 M© WD 
{dash-dotted line). 
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